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Problem: 


► The accurate computation of ro-vibrational energy levels of 
polyatomic molecules using potential energy functions 

► H2O, HCCH, ... 
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Coordinates are King! 


► Normal modes 

► KE^ complicated, BC^ nasty 

► Polyspherical Orthogonal coordinates 

► KE and BC simple, PES complicated 

► Bond-Length-Bond-Angle coordinates (Polyspherical 
Non-orthogonal coordinates) 

► Quartic PES and BC simple 

► KE complicated 
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Coordinates are King! 


► Normal modes 

► KE^ complicated, BC^ nasty 

► Polyspherical Orthogonal coordinates 

► KE and BC simple, PES complicated 

► Bond-Length-Bond-Angle coordinates (Polyspherical 
Non-orthogonal coordinates) 

► Quartic PES and BC simple 

► KE complicated 

► What is optimum? 

► BLBA: 15.43 cm“^ error (1992) 

► Normal from Symmetric Jacobi: 10.79 cm“^ 

► Radau: 1.49 cm“^ 

► Optimized: 1.29 cm“^ 
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Chosing coordinates: Step 1 


► 


X matrix of nuclear Cartesians 

► T = -^y 

2 ^1 


ax? 

IQ. 
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Chosing coordinates: Step 1 
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X matrix of nuclear Cartesians 

► T = -^y 
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X matrix of internal Cartesians (ic) 

► Normal coordinates: x = LX 

► Generalized coordinates: x = XM 


Chosing coordinates: Step 1 


X matrix of nuclear Cartesians (X) 


T = V ^ V- 
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► X matrix of internal Cartesians (x) 

► Normal coordinates: x = LX 

► Generalized coordinates: x = XM 
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Chosing coordinates: Step 1 
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► 


X matrix of nuclear Cartesians (X) 


T = -^y — V- 

2 rric 
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X matrix of internal Cartesians (ic) 

► Normal coordinates: x = LX 

► Generalized coordinates: x = XM 
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Orthogonal coordinates: 
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Role of M 
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M only enters in computation 
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evaluation of potential: X = xM ^ 


and 
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evaluation of potential: X = xM ^ 
For BLBA, = ±1, 0 
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Role of M 


► 

► 

► 
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M only enters in computation of 

evaluation of potential: X = xM 
For BLBA, = ±1, 0 

det M 0, /3 = A/ c.m. and 
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-1 



and 


. A/. 


V □ f 


<1 


t- ■« 




Chosing coordinates: Step 2 
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Chosing coordinates: Step 2 
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T =??? 


Chosing coordinates: Step 2 
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r2 cos ax 
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T =??? 




► Bramley and Handy tetra-atomic: T has 75 terms . 


Chosing coordinates: Step 2 


► X = 
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Bramley and Handy tetra-atomic: 

T has 75 terms 


basis functions: anything with finite matrix elements 


Corrections to Born-Oppenheimer Approximation (2001) 
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Corrections to Born-Oppenheimer Approximation (2001) 


► B.O.; \/ from T = 

► B.O.; V/ from , T = - f ^ 

► order: V from < + T >, T 
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Corrections to Born-Oppenheimer Approximation (2001) 
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B.O.: \/ from ^ E„ ^ E; 

icy. 

B.O.: \/ from H^'r' = E^' , 7" = -f Ea 7^ E/ ^ 

order: V from < >, T = same 

2"^^ order: changes, so V =same+C^°^ 
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Easier way to derive T 
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Easier way to derive T 
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Easier way to derive T 
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Easier way to derive T 
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Easier way to derive T 
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KJi — g^, ..., Jx, Jy, ana 

for openshell systems ... 


What If All Orthogonal Coordinates Not Very Good? 
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► SO2 

► C3 



Radau Symmetric Stretch 
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What If All Orthogonal Coordinates Not Very Good? 


► SO2 

► C3 

► C3H- 
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How to optimize M 


► Constraints: (3 = N c.m., 
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How to optimize M 
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How to optimize M 


► Constraints: (3 = N c.m. 

► M/3/3 = 1 

► V cross terms small and 

► All Tg roughly the same 


1 

M/3W 


0 


1 

M/3^' 


small, 13 ^ 13 


How to optimize M 


► Constraints: B = N c.m., = 0 

^ f^/3N 

► M/3/3 = 1 

► V cross terms small and smal 

1 ^ 13 / 3 ' 

► All Tg roughly the same 

► Repeatedly guess and walk downhi 
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Summary 


► 1980's Handy introduces Bond-length-bond-angle coordinates 
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Summary 


► 1980's Handy introduces Bond-length-bond-angle coordinates 

► 1992 I scoff at BLBA coordinates 

► 2001 I use BLBA coordinates 

► today I embrace generalized coordinates 
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